One-loop renormalization of fermionic currents with the 

overlap- Dirac operator 



C. Alexandrou a , E. Follana", H. Panagopoulos", E. Vicari b 
a Department of Physics, University of Cyprus, P.O.Box 20537, Nicosia CY-1678, Cyprus. 
b Dipartimento di Fisica delVUniversita and I.N.F.N., Via Buonarroti 2, 1-56127 Pisa, Italy. 

(February 1, 2008) 



Abstract 



We compute the one-loop lattice renormalization of the two-quark opera- 
tors TpTtp, where V denotes the generic Dirac matrix, for the lattice formulation 
of QCD using the overlap-Dirac operator. 

We also study the renormalization of quark bilinears which are more ex- 
tended and have better chiral properties. 

Finally, we present improved estimates of these renormalization constants, 
coming from cactus resummation and from mean field perturbation theory. 
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I. INTRODUCTION 



Recent developments in lattice QCD have shown that chiral symmetry can be realized on 
the lattice without fermion doubling (for recent reviews see e.g. Refs. [Jj],0), circumventing 
the Nielsen-Ninomiya theorem ||. This has been achieved by introducing an overlap-Dirac 
operator jj| derived from the overlap formulation of chiral fermions [||. The simplest such 
example, for a massless fermion, is given by the Neuberger-Dirac operator 0] 



D N = - p 

a 



1 +x{x^xy 1 



-1/2 



X = D w p, 



(1) 
(2) 



where -Dw is the Wilson-Dirac operator (with the Wilson parameter r set to its standard 
value, r = 1) 



D 



1 r 



W 



7,. (v; + v M ) - av;v, 



(3) 



V^ip(x) = -[U (x, /j,)i/j(x + ap) - if>(x)] , (4) 

a is the lattice spacing and p is a real parameter subject to the constraint < p < 2. 
Nonperturbatively one expects — m c < p < 2, where m c < is the critical mass associated 
with the Wilson-Dirac operator. 

-Dn satisfies the Ginsparg- Wilson relation || 

l5 D + D l5 = aD l5 D, (5) 

which protects the quark masses from additive renormalization . Lattice gauge theories 
with Ginsparg- Wilson fermions have been proved to be renormalizable to all orders of pertur- 
bation theory ||. The Ginsparg- Wilson relation allows us to write, at finite lattice spacing, 
relations that are essentially equivalent to those holding in the low-energy phenomenology 
associated with chiral symmetry (see e.g. Refs. |9,l0|). It indeed implies the existence of an 



exact chiral symmetry of the lattice action under the transformation (TT| 

5 e i[)(x) = ej 5 (l - aD)if)(x), S e ^(x) = ^(x^e, (6) 

which can be extended to the flavour non-singlet caseQ. The axial anomaly then arises 
from the non-invariance of the fermion integral measure under flavour-singlet chiral trans- 



formations [11-16]. A lattice formulation of QCD satisfying the Ginsparg- Wilson relation 



1 Actually Ginsparg- Wilson fermions are invariant under a more general transformation 

5 e ^{x) = e 75 (l - vaD)if)(x), 8 e $(x) = $(x) [1 - (1 - v)aD] 75 e, (7) 

where v can be chosen arbitrarily. All these symmetries are essentially equivalent, leading to 
equivalent Ward identities [|l(|. 
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overcomes the complications of the standard approach (e.g. Wilson fermions), where chiral 
symmetry is violated at the scale of the lattice spacing. The chiral symmetry ensures that 
the hadron masses are free of 0(a) discretization errors. Indeed their leading scaling correc- 
tions are 0(a 2 ). The important point is that lattice Dirac operators satisfying Eq. @ are 
not affected by the Nielsen- Ninomiya theorem ||, thus they need not suffer from fermion 
doubling. 

Indeed, avoids fermion doubling. The locality properties of in the presence of a 
gauge field are not obvious. is not strictly local, but locality should be recovered in a 
more general sense, i.e. allowing an exponential decay of the kernel of with a rate which 
scales with the lattice spacing and not with the physical quantities [T7| . Thus, the Neuberger- 
Dirac operator seems to have all the right properties that a lattice Dirac operator should 
have in order to describe massless quarks. However, its complexity renders its numerical 
implementation very demanding. In this respect some progress has been achieved (see, e.g., 
Refs. [fT7T-|27||), and Monte Carlo simulations seem to be already feasible, at least in quenched 
approximation. 

The matrix elements of the fermionic currents are employed to predict hadronic decay 
constants, electromagnetic and weak form factors, quark masses, etc. A knowledge of their 
lattice renormalization constants is necessary to relate the matrix elements computed us- 
ing lattice simulations to the corresponding ones defined using continuous renormalization 
schemes in experimental data analysis. In this paper we compute, to one loop in perturba- 
tion theory, the renormalization constants of local bilinear quark operators ipTip, where T 
denotes the generic Dirac matrix, in the lattice formulation of QCD using the Neuberger- 
Dirac operator. We also extend our computation to quark bilinears which are more extended 
and have improved chiral properties. In view of the better chiral behaviour of the Neuberger- 
Dirac operator, we expect it to have a wider use in Monte Carlo studies of hadronic physics; 
such studies require knowledge of the above renormalization constants. 

Perturbative computations are much more cumbersome than in the case of Wilson 
fermions, due to the more complicated structure of the Neuberger-Dirac operator p8|-|30|. 
One-loop calculations are already quite complicated, and require the use of symbolic ma- 
nipulation packages. We have developed such a package in Mathematica. An example of 
one-loop calculation using the Neuberger-Dirac operator is reported in Ref. ||29|| , where we 
computed the ratio Al/Aj^ between the A-parameters of the lattice formulation and of the 
MS renormalization scheme. 

Nonperturbative methods to estimate the renormalization constants, such as those of 
Refs. H3T H33H would in general be preferable to approximations based on perturbative cal- 
culations, due to their better controlled systematic errors (0(a) against 0((?q)). However, 
perturbative estimates may still be quite useful. They indeed provide important consistency 
checks. Further, in those cases where nonperturbative methods are very costly to implement, 
as in the case of the Neuberger formulation of the Dirac operator, perturbative methods may 
remain the only source of quantitative information. 

One could, of course, define exactly conserved vector and axial currents having Z4 = 
Zy = 1, following Noether's procedure, (i.e. performing variations of the action with respect 
to the non-singlet chiral transformations of Eq. (|6|)). However, these Noether currents are 
rather complicated extended objects, which might be cumbersome to use in actual simula- 
tions. Their expressions for the Neuberger-Dirac operator can be found in Ref. . 
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The paper is organized as follows: In Sec. we report the one-loop calculations necessary 
for the evaluation of the lattice renormalizations of the two-quark operators. In Sec. [II we 
improve the one-loop estimates by performing a resummation to all orders of a certain 
class of gauge invariant diagrams, dubbed "cactus" diagrams |34H36j. We also compare 



them with the results of other improvement recipes [p7|,|3§ 
improved perturbation theory. 



such as the so-called tadpole 



II. ONE-LOOP LATTICE RENORMALIZATION OF THE TWO-QUARK 

OPERATORS 

A. Formulation of the problem 

The lattice regularization of QCD we consider is described by the action 
i N f 

Sl = - £ Tr t 1 - ^Wl+EE^W^^^id/). (8) 

9o x,fi,u i=\ x,y 

where U^ u (x) is the usual product of SU(N) link variables U^{x) along the perimeter of a 
plaquette originating at x in the positive \x-v directions, and Nf is the number of massless 
quark flavours. 

The observables we study are bilinear quark operators of the form 

% = ^(x)r^(x), (9) 
where Tj denotes generic Dirac matrices, i.e. 

1, 7s, 7m> 7m7s, <V75- ( 10 ) 

Specific bilinear operators are denoted according to their Lorentz group transformations: 
S(x) = $(x)ij){x), P(x) = ^(z) 7 5 V> (a?), V^x) = ^{x^^x), A^{x) = ip(x)-f^ 5 ip(x), and 
Tn V (x) = ip(x)a fJiU 'y5'ilj(x). In the above definitions, flavour indices are left unspecified. In- 
deed, this is sufficient for our purpose, since one-loop calculations are unaffected. 

The lattice renormalization constants Zq. of the operators Oi are defined once a renor- 
malization scheme is given: we will consider the MS renormalization scheme. In order to 
simplify our calculations we work with zero quark mass. This is justified by the fact that 
renormalization constants in the MS scheme are independent of the fermionic mass. 

We also introduce the following bilinear operators 

S' = <$ (l - ±aD N ) tfj, P' ee ^75 (l - laD N ) ^ (11) 

VI ee ^ (l - |a£> N ) A 4 = ^7s (l - |«^n) i/>, (12) 

which are extended, and local only in the more general sense holding for D^. They are inter- 
esting since under the lattice chiral symmetry of Eq. (|6|) they transform as the corresponding 
operators in the continuum. Indeed one can easily check that 

5 e S' = 2P'e, S e P' = 2S'e, (13) 
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and under the non-singlet transformations 



S e a^ = e a T^ 5 (l - aD)^, 5^ = ^ 5 T»e a , (14) 

the vector and axial currents V and A' transform as 

S eb V' a = if abc e b A' c , 5 eb A' a = if abc e b V' c . (15) 

Using the corresponding Ward identities one can show that 

Z$* = Zpi, Zyi = Z A' . (16) 

We mention that, as pointed out in Refs. P,|TU|], the spontaneous breaking of the non-singlet 
chiral symmetry, Eq. (§), is related, for finite lattice spacing, to the vacuum expectation 
value (condensate) of the operator S', i.e. it occurs if (S') ^ 0. 

One can also prove that the renormalization constants of 0\ = 01^(1 — ^aDjq)ip coincide 
with those of the corresponding operators Oj = tpTiip. This can be seen from the path 
integral representation for Zq^. Upon fermionic integration, the factor of in the definition 
of 0[ will cancel against a fermion propagator; there will remain another propagator which 
will give a vanishing contribution to Zo i; since it cannot develop any 0(1/ a) singularities 
(by virtue of chiral invariance, which forbids additive mass renormalization). Thus, using 
the relations Eq. (0), we also have that 

Z s = Z P , Z v = Z A . (17) 

Note that if one introduces the bare fermionic mass by writing the Dirac operator as in 



Ref. [39 



(l - |am ) -Dn + mo, (18) 

then the renormalization of the fermionic mass is related to that of the operator 
■0^1 — \aD^ tp, denoted Z S} by the relation 

Z m = Z s \ (19) 

It is worth mentioning that the chiral symmetry of the Neuberger-Dirac operator ensures 
the absence of 0(a) discretization errors only for the spectrum of the theory. In order 
to achieve the same property for generic matrix elements of local operators, one should 
employ improved operators |4(|. This issue is discussed in Ref. |^1]] for bilinear operators 
constructed with Ginsparg- Wilson fermions; there it is shown that, for massless quarks, such 
an improvement can be achieved using the operators 

O'l = $ (l - \aD^) r, (l - \aD^) j>. (20) 

Proceeding as before, one can prove that the renormalization constants of the operators 0'( 
are the same as those of the corresponding Oj. 

The standard perturbative computation of the lattice renormalizations Zq. of the opera- 
tors Oi requires the calculation of the two-quark one-particle irreducible function, T(p), and 
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the two-quark one-particle irreducible functions ToXp) with an insertion of the operators 
Oj. In the massless case, the quark-field renormalization is obtained through the relation 



(21) 



where T L (p) is the two-quark function calculated on the lattice (in the limit a — > 0), and 
r MS (p) is the MS-renormalized two-point function, which can be computed in the continuum. 
Indeed a simple calculation gives 



r MS (p) = iPM 



l + g'c F [c\n^ + b MS \+0(g*) 



where 



a 



c = 



V 



a 



(22) 



(23) 



16tt 2 ' 16tt 2 ' 
and a is the gauge parameter. The renormalizations Z . can be obtained by the equation 

rgf(p) = Z Xa^Z4a^T L 0i (p) (24) 

where T (p) and Fq S (j)) are the two-quark functions with an insertion of calculated 
respectively on the lattice and in the continuum using the MS renormalization scheme. 
Setting 



MS,L, 



one has 



B^(P) = 1 + 9 2 c F I c Ql In ^ + e y ] + 0(/), 



/r 



where (setting a = 1) 



cs,p 



4tt 2 ' 



and 



lMS _ J 



Cv,A 



h MS _ 



167T 2 ' 



v,a 16n r 



c T = 0, 



b™ s = 0. 



(25) 



(26) 



(27) 



(28) 



In the above formulae we have neglected terms proportional to p^Pvjp \ these are present 
also in the lattice expressions, and cancel out in the calculation of the renormalizations. 
On the lattice one finds 



and 



r L (p) = ip^ [l + g 2 c F (-clnoV + b L ) + 0{g A ) 



BUp) = ±+9 2 cf (-co, lna 2 p 2 + b L 0i )+0(g 4 ). 



(29) 



(30) 
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Using Eqs. (pT|) and (p4|), we obtain 



Z i> = \ + g 2 c F (c In aV + & MS - & L j + 0(<? 4 ), 
Z Gl = 1 + <?V [(co, - c) In a 2 /? + 6of - ^ - ^ + & L 1 + 0(^ 4 ) 



MS 



(31) 
(32) 



Note that Zo t are independent of the gauge parameter a. 

So, to compute Zq. we need to evaluate the constants b L and 6q. by a one-loop perturba- 
tive calculation on the lattice. We have to expand the Neuberger-Dirac operator in powers 
of g , and use the resulting vertices to construct the diagrams related to the one-particle 
irreducible functions T L (p) and Tq.(p). In the appendix we list the relevant formulae for our 
one-loop calculations. 



B. Results and discussion 



Two diagrams contribute to T L (p) , shown in Figure 1 ; for Bq . (p) , there is only one dia- 
gram to one-loop, shown in Figure 2. Given that the 4-point vertex contains a part with an 
internal momentum (k in Eq. ( |A9| ) ) , the corresponding part of the second diagram in Fig. 
1 actually has the same connectivity as the first diagram. 



(1) 




(2) 




FIGURES 1, 2. One- loop diagrams contributing to the functions T L (p) (Fig. 1), and Bq.(p) (Fig. 
2). Wavy (solid) lines represent gluons (fermions); a cross is an insertion of the operator Oi . 
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The algebra involving lattice quantities was performed using a symbolic manipulation 
package which we have developed in Mathematica. For the purposes of the present work, 
this package was augmented to include the propagator and vertices of the overlap action. 

We express our results in the form of Eqs. (|29|) , flSCf ). The parameters 6 L , 6q. depend on 
p, but not on N or Nf. 

To extract the p-dependence, we first isolate the divergent terms; these are responsible 
for the logarithms. There are only a few such terms, and in the pure gluonic case their 
values are well known. We can use these values also in diagrams with fermions, applying 
successive subtractions of the type: 



Q 




(33) 



where q 2 is the inverse fermionic propagator. All remaining terms now contain no diver- 
gences, and can be evaluated by Taylor expansion in ap. 

At this stage, one is left with expressions which no longer contain p and must be nu- 
merically integrated over the loop momentum. Given the complicated form of the overlap 
vertices, these expressions turn out to be quite lengthy, containing a few hundred terms in 
the cases at hand. 

The integration is done in momentum space over finite lattices; an extrapolation to in- 



finite size is then performed, in the manner of Ref. ||42|| . We evaluated the integrals for a 
range of values of the parameter p, as presented in Table 1. For all values of p that we 
quote, lattice sizes L < 128 are sufficient to yield answers to at least 7 significant digits 
(the uncertainty coming from a systematic error in the extrapolation, which can be esti- 
mated quite accurately). As the endpoints of the perturbative domain of p are approached 
(p — > 0, p — > 2), some of the quantities we calculate require increasingly larger lattices, for 
similar accuracy; this is, of course, a reflection of the divergences in the propagator at these 
endpoints. Figure 3 shows the dependence of our results on p, for the whole range < p < 2. 

A number of consistency checks, some of which are rather nontrivial, may be performed 
on our results. In particular: 

• The logarithmic coefficients must equal those of the continuum. 

• Terms proportional to p^Pu/p 2 , which appear in Z v , Za, should match those of the 
continuum. 



• No O(p ) terms must appear in T L (p), i. e. no additive mass renormalization, as 
required by chiral symmetry. 

Our results fulfill all of the above requirements. The last one serves also to verify our 
estimates of the systematic errors coming from the extrapolation. 
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FIGURE 3. The coefficients b (p), bQ.(p), as a function of the parameter p . 



For quick reference, we write the one-loop values of Z^(a = 1), Zq 1 at p = 1, as follows: 



= 1 + g 2 c F [(\na 2 fi 2 )/167i 2 - 0.231966] 
Z S>P = 1 + g 2 c F [3 (lnaV)/167T 2 + 0.204977; 



j a,v 



1 + g 2 c F [ 0.198206; 



Z T = 1 + g 2 c F [ -(lnaV)/16vr 2 + 0.204392 



(34) 
(35) 
(36) 
(37) 



III. IMPROVED PERTURBATION THEORY 

In order to improve the estimates coming from lattice perturbation theory, one may 
perform a resummation to all orders of the so-called "cactus" diagrams |34H36l. Briefly 
stated, these are gauge-invariant tadpole diagrams which become disconnected if any one 
of their vertices is removed. The original motivation of this procedure is the well known 
observation of "tadpole dominance" in lattice perturbation theory. In the following we refer 



to Ref. 1 34] for definitions and analytical results. 
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Since the contribution of standard tadpole diagrams is not gauge invariant, the class 
of gauge invariant diagrams we are considering needs further specification. By the Baker- 
Campbell-Hausdorff (BCH) formula, the product of link variables along the perimeter of a 
plaquette can be written as 

= exp \ig (A x>fJt + A x+tM>u - A x+V ^ - A XjU ) + 

= exp {ig F£l + ig 2 F%l + O^) } (38) 

The diagrams that we propose to resum to all orders are the cactus diagrams made of 
vertices containing F^L, ■ Terms of this type come from the pure gluon part of the lattice 
action. These diagrams dress the transverse gluon propagator Pa leading to an improved 
propagator P4 , which is a multiple of the bare transverse one: 

PT = J^T-y (39) 

i - w{go) 

where the factor u>((?o) will depend on g and N, but not on the momentum. The function 



w(g ) can be extracted by an appropriate algebraic equation that has been derived in Ref. p4 
and that can be easily solved numerically; for SU(3), w(g ) satisfies: 



ue- u/3 



M 2 /3 _ Au + 8 ] = 2g 2 , u(g ) = 9 % ■ (40) 



The vertices coming from the gluon part of the action, Eq. (||), get also dressed using a 
procedure similar to the one leading to Eq. fl39|) [34]. Vertices coming from the Neuberger- 



Dirac operator stay unchanged, since their definition contains no plaquettes on which to 
apply the linear BCH formula. 

One can apply the resummation of cactus diagrams to the calculation of the renormal- 
ization of lattice operators. Approximate expressions are obtained by dressing the cor- 
responding one-loop calculations. Applied to a number of cases of interest |34]j35| , this 



procedure yields remarkable improvements when compared with the available nonperturba- 
tive estimates. As regards numerical comparison with other improvement schemes, such as 
boosted perturbation theory |37|j3^| , cactus resummation fares equally well on all the cases 



studied |36]]. It is worth mentioning in passing that cactus resummation also affords us a 
systematic means of improving perturbation theory, by successively dressing higher loop 
diagrams. 

Let us consider the renormalizations Zy, Za of isovector fermionic currents, which are 
finite functions of the bare coupling go- At one-loop order we have 

Z VA = 1 + gfiz v ,A + ■■■, (41) 

where the constants zy and za have been calculated in Sec.||. The cactus dressing of the 
above one-loop expressions can be simply obtained by using the dressed transverse gluon 
propagator, Eq. (|3~9l). We thus obtain the following approximate expressions 
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Let us apply the above formula to the lattice SU(3) gauge theory, for which Zy >A = 0.26427, 
at go = 1 which is a typical value for Monte Carlo simulation. Since |5jJ], 



1 - w ( go = 1) = 0.749775, 



(43) 



we find Z VtA (p = 1) ~ 1.35 

Other recipes of improvement have been proposed in the literature (see e.g. | 3*5fl , and |43 
for a review of them) that essentially consist in a better choice of the expansion parame- 
ter. Among them we mention the so-called tadpole improvement [38[ (MFI) motivated by 



mean-field arguments, in which one scales the link variable with u (go) = (jj^ 1 U x ^ v ) 1 ^ as 
measured in the Monte Carlo simulation. Accordingly one rescales the coupling constant: 
9o ~^ 9mf = 9o/ u o- Thus, one obtains a mean-field improved expansion 



'V,A = U 



1 + 9 



mf 



z v ,A + — )+0(g 



4 1 

mfJ 



(44) 



For example, for 577(3) in the quenched approximation and at g% = 1 one finds uq ~ 0.878 
and g^ f ~ 1.68. Putting these number in Eq. ( f44|) we obtain Zy,A — 1-39, which is in 
reasonable agreement with the estimate coming using the cactus resummation. 
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APPENDIX A: 

In order to perform the lattice perturbative calculation we must formally expand in 
powers of go- We list here the relevant expressions for the propagator and vertices, following 



Ref. [28 



Let us first write down the weak coupling expansion of the Wilson-Dirac operator Df. 
This will be useful for constructing the relevant vertices of . We write 

X{q,p) = D w (q,p) --p = Xo(p)(27rY5 A (q-p)+X 1 (q,p)+X 2 (q,p) + 0(g 3 ), (Al) 

(X 



where 



i 1 1 

x o(p) — - ^2ifi s ^ na p^ + _ XX x — cos a P») — p> ( A2 ) 



a „ a ~r a 



Xi(q,p) = g Q J d^iq-p-^A^V^ip + k^), (A3) 
Vi,p(q) = ^cosag^ + sinag^, 
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MQ,P) = f / d4 (2ny 2 6{q - p - kl ~ ^MMW^ip + h/2 + k 2 /2), (A4) 
V2,»{q) = -?7m° sin a % + ° cos a %- 
The Fourier transform of the Neuberger-Dirac operator takes the form 
1 



■D N (q,p) = D (p)(27iY5\q-p) + Z(q,p) 



D (p) is the tree level inverse propagator: 

D \P) 



2[cu(p) + b(p)} 2' 



where 



w (p) = - f 1] sin2 + E(! - cosap M ) - pf j 

8 \ (J fl / 



1/2 



= -ZK 1 - C0SCI Pm) - -P- 



The function can be expanded in powers of g as 



Lu(p)+u(q) 
1 



+ 



u(p)+u(q) 
d 4 k 



*i(<Z,J>) - 
X 2 (q,p) 



u{p)u{q) 
1 



/ 



u)(p)u(q) 
1 



X (p)Xl(p,q)X (q) 
X (p)Xl(p,q)X (q) 



x 



(A5) 



(A6) 



(A7) 
(A8) 



(2tt) 4 uj( P ) + w(p) + uj(k) u(q) + 
-X Q (p)Xl(p, k)Xi(k, q) - X 1 [p, QXtWX^k, q) - X,(p, k)X\(k, q)X (q) 
u(p) + u(q) + uj(k) 



u(p)u(q)u(h) 



-X (p)Xl(p, k)X (k)Xl(k, q)X (q) 



+ ... 



(A9) 



From Yl(q,p) one can read off the vertices necessary for the one-loop calculations presented 
in this paper. 
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TABLE I. Coefficients b L (p), b^. (p). 
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b L 


b L 


°V A 


b^ 


0.2 


1.5236150 


0.0450336 


0.0399790 


0.0340724 


0.3 


0.9841497 


0.0486865 


0.0399888 


0.0328679 


0.4 


0.7155175 


0.0517536 


0.0399997 


0.0318600 


0.5 


0.5550679 


0.0544426 


0.0400116 


0.0309796 


0.6 


0.4486267 


0.0568654 


0.0400248 


0.0301895 


0.7 


0.3729964 


0.0590910 


0.0400393 


0.0294670 


0.8 


0.3165894 


0.0611655 


0.0400554 


0.0287970 


0.9 


0.2729750 


0.0631217 


0.0400733 


0.0281687 


1.0 


0.2382987 


0.0649842 


0.0400931 


0.0275744 


1.1 


0.2101107 


0.0667722 


0.0401153 


0.0270079 


1.2 


0.1867796 


0.0685012 


0.0401400 


0.0264646 


1.3 


0.1671775 


0.0701842 


0.0401679 


0.0259407 


1.4 


0.1504998 


0.0718326 


0.0401993 


0.0254331 


1.5 


0.1361573 


0.0734568 


0.0402349 


0.0249392 


1.6 


0.1237088 


0.0750665 


0.0402755 


0.0244568 


1.7 


0.1128172 


0.0766713 


0.0403222 


0.0239840 


1.8 


0.1032209 


0.0782815 


0.0403762 


0.0235194 
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